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Abstract 

The phase structure of ABJM theory with mass m deformation and non-vanishing 
Fayet-Iliopoulos (FI) parameter, (, is studied through the use of localisation on § 3 . 
The partition function of the theory then reduces to a matrix integral, which, in 
the large N limit and at large sphere radius, is exactly computed by a saddle-point 
approximation. When the couplings are analytically continued to real values, the 
phase diagram of the model becomes immensely rich, with an infinite series of third- 
order phase transitions at vanishing Fl-parameter [I]. As the FI term is introduced, 
new effects appear. For any given 0 < £ < m/2, the number of phases is finite and 
for f > m/2 the theory does not have any phase transitions at all. Finally, we argue 
that ABJM theory with physical couplings does not undergo phase transitions and 
investigate the case of?7(2)x[/(2) gauge group in detail by an explicit calculation 
of the partition function. 
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1 Introduction 


In some cases, in supersymmetric gauge theories, observables with enough supersymmetry 
may be computed exactly. For example, this applies to the partition function in J\f = 2 
four-dimensional gauge theories on S 4 through the technique of localisation |2|. This 
allows us to write observables with a sufficient amount of supersymmetry, such as the 
partition function, in terms of matrix integrals. Even though being much simpler than 
the original functional integrals, these matrix integrals still carry important information 
of the underlying held theory. 

Under favourable conditions, these integrals may sometimes be determined exactly, 
though this often amounts to going to the large N, or planar, limit. In this limit, one 
may use techniques from random matrix theory to solve these matrix integrals |3|. One 
of the most interesting applications of this approach is perhaps the test it provides of 
the conjectured gauge/gravity-duality: localisation and large A r -techniques have provided 
insights into the strong-coupling behaviour of held theories with holographic duals, thus 
allowing for direct comparisons of quantities with a non-trivial dependence of the coupling 
constant between the gauge- and gravity-side of the duality, with excellent agreement (for 
a review, see [4]). 

In recent years, the occurrence of quantum phase transitions in massive gauge theories 
in the decompactihcation limit (where the radius of the sphere is taken to infinity), has 
been found in a wide variety of supersymmetric theories. They were hrst found for 
J\f = 2 four-dimensional gauge theories with massive matter |[5], whose critical properties 
were further investigated in 

BE 3 ® 


6||n|; 

and hve dimensions 1151 


and, soon after, new examples were found in three- 
A generic type of phase transition occurs when, at 


specific couplings, extra massless states appear in the spectrum and contribute to the 
saddle-point. It is a resonance phenomenon, as explained in |5j. These resonance effects 
already appear in theories with fundamental matter, but become much more complicated 
when the theory contains adjoint matter, leading to an infinite sequence of secondary 
resonances and a much richer phase structure. In four dimensions, the simplest example 
with massive adjoint matter is the M = 2* theory, which, in the decompactihcation limit, 
exhibits an infinite series of weak/strong quantum phase transitions as the coupling grows. 
Moreover, these phase transitions accumulate at strong coupling [5]-7|, raising interesting 
questions about how they manifest in the holographic dual theory. These were studied 
in detail in [8,10 , and similar phase transitions are also found to be present when the 


four-sphere is squashed |9|. Recently, it was shown that these non-trivial phases induce 
non-analytic dependence of higher-rank Wilson loops on the rank parameter 12 , which 


could in principle have a holographic counterpart in terms of phase transitions in the 
effective held theory on the dual D brane description. 

This type of phase transitions were also shown to occur for low-rank gauge groups, 
such as M - 2 SU(2) gauge theory with fundamental matter 11 , i.e. supersymmetric 


577(2) QCD with two massive flavours. In this case, the quantum critical point of these 
phase transitions is identified with the Argyres-Douglas superconformal fixed point of 
the theory. For SU(2) gauge group, the phase transition is driven by instantons (unlike 
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the large N case, where they are negligible), and, in the decompactihcation limit, the 
free energy log Z including instantons can be expressed in terms of the Seiberg-Witten 
prepotential -a connection that arises when a saddle-point exists 11 . 


Similar phase transitions occur in three-dimensional Chern-Simons theories with mas¬ 
sive matter. For example, three-dimensional Chern-Simons theory with massive funda¬ 
mental matter has been shown to undergo three distinct phases as the coupling runs from 
zero to infinity 13 . These phase transitions have furthermore been investigated at finite 


N, where the finite N partition function of U(N ) Chern-Simons theory with massive 


fundamental matter was computed using the theory of Mordell integrals 14 


The theory which will be in focus in this paper is the three-dimensional ABJM theory, 
which, in the massless case is an J\f - 6 superconformal theory dual to type IIA string 
theory on AdS 4 x CP 3 [16 . By using localisation techniques on the held theory side, 
the partition function can be computed exactly at large N 17,18 , which may then be 


compared to geometric analyses on the string theory side, enabling extensive tests of the 
gauge/gravity conjecture 18,19 . The partition function can be computed exactly even 


in the mass-deformed case as the path integral localises to a matrix model on § 3 17,20 


Recently, it was shown that a version of ABJM theory obtained by analytic continuation in 
the couplings exhibits phase transitions [l]. This theory contains bi-fundamental matter 
and the resulting phase structure arising in the decompactihcation limit resembles the case 
of four-dimensional M - 2* theory, with an infinite number of third-order phase transitions 
accumulating at strong coupling. These models corresponds to computing the partition 
function in the region of parameter space where the couplings are real (and the Chern- 
Simons levels are imaginary). This method, where one starts with unphysical couplings, 
was successfully implemented in the past [4,18 to compute the large N partition function 
for the ABJM theory with no mass deformation. A direct, analytic calculation of the 
large N partition function in ABJM theory with physical Chern-Simons levels is more 
complicated and has not been carried out so far. The reason is that eigenvalues appear to 
be distributed in different cuts in the complex plane with non-homogeneous N dependence 
for real and imaginary parts (see e.g. a discussion in [19]). 

In [l], the Fayet-Iliopoulos (FI) parameter was set to zero. The aim of this work 
is to incorporate this parameter, which, as we shall see, significantly enriches the phase 
structure of the theory. As the Fl-parameter approaches zero, the phase structure found 
herein agrees with previous results, and furthermore serves to clarify some of the pecu¬ 
liar behaviour previously seen. However for non-vanishing Fl-parameter, the situation is, 
as mentioned, significantly different, and the solutions to the saddle-point equations are 
divided into cases depending on the precise relations between the mass-deformation pa¬ 
rameter m and the Fl-parameter, (. The most dramatic effect is that the number of phase 
transitions undergone as the coupling is increased are now finite, and for ( > m/2 there is 
only a trivial phase with constant eigenvalue density in some finite region of support. 

The structure of this work is as follows: In section [2j we give the matrix model 
representation of the partition function for AB J (M) theory with mass and FI deformation 
parameters, and discuss some basic properties. In section [3j the analytical continuation is 
introduced. In subsection [3Tl we specialise to the situation where the couplings Ai, A 2 for 
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the two (analytically continued) gauge groups are equal and real, and in sections [I] and [5j 
we determine the eigenvalue densities for the theory in the large N limit. As mentioned, 
the precise form of these depend on the relations between m and and these different 
cases are considered in sections [4] and [5] In section [6j we provide a summary of our results 
for the cases of equal and real couplings, and also discuss the case of Ai + A 2 ; A i)2 > 0, 
for the two different schemes of analytic continuations. In section [7| we present a general 
argument showing that physical ABJM theory does not have phase transitions. The 
absence of phase transitions is illustrated by considering the case with U( 2) x 17(2) gauge 
group, where the partition function can be explicitly computed both for the analytically 
continued model and for the model with physical couplings. The analytically continued 
model exhibits analogous phase transitions as in the large N case. However, the physical 
ABJM theory does not undergo any phase transition. One key feature that seems to make 
physical ABJM theory to be non-generic is the fact that the Chern-Simons levels for the 
gauge groups are equal and opposite. Finally, in section [8j we end with some concluding 
remarks. Details of the calculation of section |5.4| are given in an appendix. 


2 Deformed ABJ(M) theory 


ABJM theory is a three-dimensional superconformal theory with maximal supersymmetry, 
gauge group Uk(N) x U-k(N), and matter in the bifundamental representation, k denotes 
the Chern-Simons level of the two gauge groups respectively. In ABJ theory 21 , the 


situation is generalised so that the two gauge groups are allowed to have diffterent ranks. 
Allowing for this small generalisation was shown to be useful when computing the partition 
function of the massless theory by analytic continuation 
deformed case [I]. 

In the massless case, these theories are well-studied and their path integrals on § 3 are 


18 , as well as in the mass- 


known to localise onto constant held configurations 17,22 . However, these results do not 


rely on the theory being conformal, nor maximally supersymmetric, and thus they may 
be used to examine the deformed ABJ theory as well. 

In previous work, the decompactihcation limit of the mass-deformed ABJ theory was 
considered [l], and quantum weak/strong phase transitions were shown to be present in 
two analytically continued versions of the theory. 

However, there is another way of introducing a scale to the problem, other than in¬ 
troducing a mass: by introducing a Fayet-Illiopoulos deformation. This will be present 
in the most general form of supersymmetric deformation, and it is this setup which will 
be considered herein, with both non-vanishing mass and Fl-parameter. The partition 
function on § 3 takes the form of an eigenvalue integral in the large 77-limit 17,20 , and, 
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with the normalisation of 18 , may be written as: 


] r N 1 rjn. Jv 2 fill 

Z ABJ m( 2C, m; fc) = f 


( 1 ) 


n sinh 

i<j 


2 Pi~Pj 
2 


[I sinh 


a<b 


2 ^g~^b 


-EC Ew+2>. 


n cosh cosh 




where p t ,zy, represent the eigenvalues of the auxiliary fields from the vector multiplets of 
the two gauge groups, m denotes the mass-deformation and ( represents the Fl-parameter, 
which will be taken to be real. 

By shifting the integration variables, x = p, - <f,y = v + (, we may then move the FI- 
dependence from the exponent to the denominator, and what we are left with may be 
thought of as the mass-deformed case but with two different masses, mi and m 2 . In these 
variables, the partition function becomes: 




C 2 (N?-N, ) 

e 2 ° rr d Va 

Ni'.Nf. J 2,n 2 tt 

n Sinh 2 ^ n Sinh 2 » 

i<j ~ a<b 

X n cosh cosh Xi-Va-m2 - e 



( 2 ) 


where g = -p represents the coupling, and mi,m 2 relates to m and ( as: 

mi - m + 2( and m 2 = m - 2( . 


(3) 


Consider the partition function (|Tj) with Ni - N 2 = N. It can be written in another 
form, which is useful to exhibit some symmetries. This is done by following 20 , slightly 
generalising their derivation for k - 1 to arbitrary k. By using the identity, 


n,;<j sinhfo - Xj) sinh (y t - yf) 




(4) 


riij cosh(a: i - yf) ^ V cosh (a;* - y p ^) ’ 

where p runs over all permutations of {1,..., N}, the partition function may be written 


as: 


-^abjm(2C,^; k) - ^ 


p,p 


(-1 )p + p' f d N p d N v 
N\ 2 

(-1 y 


i 


n N\ J 


(2 tt) n (2tt) n 
d N n d N u 


n 

i 

n 


e -iC (/A+^+fl/y 2 -^ 2 ) 
cosh cosh w ~ v 2 (<) ~ m 

e -&«w+«'0+g0M) 


(2n) N (2n) N V cosh cosh 


We now make use of the Fourier transform 

0 ITfl 


I 


dr- 


cosh TIT 


cosh | 


(5) 


( 6 ) 
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for all hyperbolic cosines in the denominator, introducing new integration variables r t , t[. 
The integrals over p*, zy then become Gaussian and can be computed explicitly, after 
which one finds 


Zabjm( 2C, m; k) = Y, ppp J d N rd N r' 


1 Zi T i( T i-T p (i))+iT, i (T.m 1 -T i m2) 

Yli cosh(7TTj) cosh(7rr/) 


(7) 


Using again the Fourier transform ([6]) and computing the integral over r/, we hnally obtain 
(upon rescaling r* -> /cr*) 


Z ABJM (2(,7n-,k) = XX _1 ) P pj f d ' Nr 


e -ikm 2 Ei n 


Fli COSh(/c7TTj) COsll(7r(Xj - Tp(j)) - 1J f) 


( 8 ) 


We stress that the derivation above only holds when the Chern-Simons levels of the 
two gauge groups U(N) kl x U(N) k2 are opposite, k 2 = -ki, which is the case in ABJM 
theory (for k 2 ± ~k \, terms t/ 2 ,t 4 2 remain in the exponent, leading to more complicated 
expressions). 

The partition function (|TJ) has the obvious symmetry ( -* -((, under which m\ •«-> m 2 . 
However, the partition function written in the form (|8]) makes manifest another symmetry, 

m 2 -m 2 , (9) 


arising after the sum over permutations. Under this symmetry, the FI- and mass-deformations 
are exchanged. I 11 other words, the deformed ABJM partition function (|T|) with N\ - N 2 
enjoys the property 

Zabjm( 2C, m; k) = Z AB jm("p2C; k) . (10) 

In particular, a Fl-deformation on the massless theory ( = m/2 is equivalent to a mass- 
deformation m in the theory with vanishing Fl-parameter, 


Z ABm (m, 0; k) = Z A BJM(0,m; k) . 


( 11 ) 


3 Analytic continuation and saddle-point equations 

Following [I], we first assume independent Chern-Simons levels k\ and k 2 for the two 
gauge groups U(Ni) and U(N 2 ), and introduce two different couplings, 

27ri 27ii 

9 1 = — ) 9 2 = — , 

k\ k 2 

and equivalently, the two’t Hooft couplings 

Ai = N\(ji , A2 = N 2 g 2 , 


( 12 ) 

(13) 


for the different gauge groups. 
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Our starting point will be the representation (J2]) for the partition function. In [I], 
two analytic continuations of this model to arbitrary Ai, A 2 were considered: one where 
we set &2 = -k\ - k and leave N \, N 2 arbitrary, which may be thought of as an analytic 
continuation in the gauge group rank. The other one may be thought of as a continuation 
in the Chern-Simons level instead (and holds the rank of the two gauge groups equal). We 
will mainly use the first analytic continuation (which is the one used in 4,[l8j| for ABJM 


theory at large N). The second analytic continuation will be discussed in sections 3.2 


6.2.2 and later in section 7.2 for the special case where the gauge group is U( 2) x 1/(2). 

It is important to note that, after analytic continuation to two independent couplings 
Ai, A 2 , the resulting partition function Z(2£,m; Ai,A 2 ) cannot be written in the form (|8|) , 
except in the special case A 2 = -Aj. I 11 particular, for generic Ai, A 2 , the partition function 


does not satisfy the symmetry (10). 


The ABJM partition function (J2| with integer k is given by a convergent integral, 
therefore in principle one does not need to resort to analytic continuation to define it. 
However, for integer k (thus imaginary A 2 = — Ai), the saddle-points lie in cuts in the 
complex plane which are complicated to determine even numerically. Here we perform 
analytic continuation to real, positive couplings Ai,A 2 because it is in this case that the 
saddle-point equations can be solved explicitly in terms of closed formulas. 

With the analytic continuation in the gauge group rank, we set = -k± = k, and the 
saddle-point equations of (J2]) take the form: 


Ai Xi-Xj A 2 vW , , x t - y a + mi Xi-y a -m 2 

x, =— > coth-+-> tanh -+ tanh- 

A '1 fy, 2 2 JV 2 M 2 2 


A 2 ^ y a - y b . A, 

A r 2 bia 


/v 2 y-' , 1 ya yb m / , . Ua x i m l , , Ua + m 2 

y n =— > coth-+-> tanh-+ tanh- 

ya t\j_ z_j 0 2 N, V V 2 2 


) (M) 

)■ 


The ABJM theory is recovered by analytic continuation Ai e^X, X 2 -> e^X, where ip 
goes from 0 to 7t/2. 

Phase transitions typically occur in the decompactihcation limit, where the radius R 
of the three-sphere (set to unity in previous formulas) is sent to infinity. The dependence 
on the radius can be restored by rescaling m -»• mR , ( -»• £ R , Xj -»■ XiR, yi ->• yiR. For 
the coupling, we take the same scaling used in |JT, 13 , where X/R is fixed as R -»• 00 . This 
particular decompactihcation limit turns out to be self-consistent and the dependence 
on R completely cancels from the saddle-point equations. As R -»• 00 , the hyperbolic 
functions are replaced by sign functions. Furthermore, in the large N limit, the eigenvalues 
Xi, y a have continuum distributions described by unit-normalised eigenvalue densities, 
p x (x)j Py(y)j an d the saddle-point equations take the form 

x =Ai dx'p x (x')sign(x - x r ) + ^ J dy p y (y){s\gi\(x - y + m{) + sign(x - y - m 2 )), 

( 15 ) 

V =A 2 J c dy' p y (y')sign(y -y') + y £ dx p x (x)(sign(y -x-mi) + sign (y -x + m 2 )), 
















where C x and C y denote the intervals on which p x and p y are supported respectively. 

It will be convenient to introduce another change of variables, namely the rescaling 
y -* -y. By also defining p y (y) = p y (-y ), such a rescaling leads to equations which are 
symmetric under the exchange of x ++ y (and as such also p x +->■ p y and the integration 
regimes) together with Ai <-»• A 2 . These equations may explicitly be written down as: 


x = 


Ai J dx'p x (x')sign(x - x') + y J dy p y (y)(sign(x + y + mi) + sign(x + y - m 2 )), 

(16) 

y =A 2 dy' p y (y')siga(y - y') + y dx p x (x)(sign(y + x + m l ) + sign (y + x- m 2 )). 

Take C x = [~A,B] and similarly C y = \-C,D], where {A,B ,C,D} € R. Differentiating (16) 
with respect to x,y respectively gives us: 


P x (-) ~ mi ) + Pv(~ z + m 2 )), 


(17) 


Pv(z) - Av{p^~ z _ mi ) + px (~ z + m ^)- 


2A 2 2A 2 

Solving these coupled functional equations is very complicated in the general case. For 
simplicity, in most of our discussion, the situation of equal, real couplings will be consid¬ 
ered. In section 6.2, we will also treat the case of generic Ai ) 2 > 0 and show that it exhibits 
the same qualitative features. 


3.1 The case of equal, real couplings 


By symmetry, it is clear from the expression of (17) that, for Ai = A 2 = A, the system 
admits a solution with two equal densities p x ,p y ■ The problem thus reduces to finding 
the solution to one single equation for a density p(z): 


p(z) =— - p(- 

’ 2A 2 ry 


■mi) - -p(-z + m 2 ), 


(18) 


where p(z) is supported on some interval \-A,B] along the real axis. 

However, unlike the case previously considered where the Fl-parameter vanishes [I], 
there is no reflection symmetry around the origin of these equations, and we cannot 
assume p(-z) = p(z). This complicates the situation compared to the case studied in [I]. 
With no loss of generality one can take ( > 0 and m > 0. Then m\ will always be greater 
than zero whereas m 2 € [- 00 , m\. 

It is clear that the solution to this equation will behave qualitatively different de¬ 
pending on the sign of m 2 , and we may thus divide our investigation into two separate 
cases: 


• rrii > m 2 , m 2 < 0 corresponding to £ > y 

• mi > 0, m 2 > 0, corresponding to ( < ■ 

These will be considered in sections [4] and [5] respectively. 
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3.2 Phase transitions in Chern-Simons theory with massive ad¬ 
joint matter 

The saddle-point equations for the second analytic continuation may be expressed in terms 


of gi,g 2 an d N, and will in this notation differ from the equations (15) of the first analytic 


continuation by some signs. Defining aq = Ngi, a 2 = Ng 2 , these equations may be written 


as: 


= J c dx'p x (x')sign(x-x') - ^ J c dyp y {y)[s\gn(x -y + mi) + sign (x-y -m 2 )^j (19) 
J c dy'p y (y')sign(y - y') f c dxp x (x){siga(y - x - mi) + sign(p - x + m 2 ))■ 

e^oq, a 2 -»• e _ *^aq, with <p 


aq 

y_ 

ol 2 


ABJM theory is recovered by analytically continuing a\ 
varying between 0 and 7r/2 . 

Consider the particular case aq = a 2 = a. We are led to a single equation: 

= 2a + \ pt '~ X \ p (~ x * ' 


( 20 ) 


If we further assume that mi = m 2 = m (i.e. ( = 0), we have reflection symmetry, and the 
equation becomes 


p( x ) ~m) + i p{x + m) . 


( 21 ) 


The solutions to this equation were studied in section 4 of jl]. 

It is interesting to compare this equation with the saddle-point equation that arises in 
a closely related system, J\f - 3 supersymmetric U(N ) Chern-Simons gauge theory with 
two massive adjoint nmltiplets. This is a precise three-dimensional analog of the J\f - 2* 
theory whose critical properties were studied in (5j-f8,10 . The partition function can be 


constructed with the general rules given in 17 , whereupon one obtains: 

N dp-i 


z -wJ ntn 


sinh 2 ^ 




iKj: cosh cosh 


( 22 ) 


At large N (and in the decompactification limit), this partition function can be computed 
by solving the same saddle-point equation (21), with a = gN. Therefore, the results of [I] 
equally apply to this case, and they can be summarized as follows: The theory has an 
infinite sequence of phase transitions, where in each phase the eigenvalue density is given 
by a piecewise constant function. The number of discontinuities increases whenever the 
coupling crosses critical values taking the theory into a new phase. In the strong coupling 
limit, a » 1, the phase transitions accumulate and the equation becomes differential, 
-m 2 p" = 1/a and p(x) approaches the asymptotic form jl] 


Poo 0*0 = 


2 grri- 


-(pg-/i 2 ), p 0 = 


/3 qm 2 \ 

— 


(23) 
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In the infinite a limit, this smooth, parabolic asymptotic density arises as the envelope of 
(discontinuous) piecewise constant densities. In the case of J\f - 2* theory, at strong cou¬ 
pling, the eigenvalue density also reaches a smooth asymptotic form after going through 
an infinite number of phase transitions, each phase described by a discontinuous density. 
In this case, the asymptotic density has the Wigner’s semicircle shape, a property that 
has been matched with the holographic prediction 


23,24 


In conclusion, J\f = 3 supersymmetric Chern-Simons gauge theory with two massive 
adjoint multiplets has large N phase transitions, with a behaviour that shares similar 
features as its four-dimensional relative, J\f = 2* theory. 


4 Equal real couplings with £ > y 


Let us now move back to the problem of ABJM theory analytically continued in the gauge 
group rank in the special case of equal and real couplings for the two gauge groups, as de¬ 
scribed by equation (18). Furthermore, let us specialise to the case of ( > y, corresponding 


to mi > 0, whereas m 2 < 0. 

Let us first consider two limiting cases, namely ( -*■ 0 and £ -*■ 00 . In the first one 
of these, the situation reduces to the one considered in jT], and equation (18) becomes 
symmetrical under reflection through the origin, implying A = B. As ( goes to infinity, 
equation (18) reduces to p(z) = ^t, since the two shifted terms will vanish as their argu¬ 


ments will lie outside the region of support for p. Knowing the eigenvalue density, the two 


interval endpoints may be obtained from the integral equation (16) with x = B, together 
with the normalisation condition. 


From these, we find A - 0, B - 2A for ( -»• 00 . It is 
then natural to assume that B > A at intermediate values (this is also confirmed by the 
numerical solution). 

With this assumption, together with the condition m\ > -m 2 > 0, the integral equation 


(16) with x - B gives 


B = 2\ 




(24) 


It can then be shown that the other interval endpoint will always lie in the origin, i.e. 
A = 0. 

This gives us an eigenvalue density as 

p(A) = ^ Vz € [0,2A] . (25) 

Therefore, there are no phase transitions in this regime where ( > y. 

I 11 conclusion, turning on a Fl-parameter ( > ™ implies a theory free from phase 
transitions. In particular, the masslesss theory with only FI deformation does not have 
phase transitions. 
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5 Equal real couplings with ( < y 


As the Fl-parameter decreases below m/2, however, the situation becomes more compli¬ 
cated: m 2 changes sign, becoming strictly positive. When ( - 0, one has m i = m 2 , which 
is precisely the case discussed in [[]. As previously mentioned, the saddle-point equation 
(18) then has reflection symmetry, which gives us that A - B, i.e., the eigenvalue density 


is supported on the interval [ -B,B ]. On the other hand, as shown above, for all £ > -?, 
the leftmost interval endpoint lies at the origin. As ( increases from 0 to the leftmost 
interval endpoint, - A , thus moves from -B to the origin. 

Precisely how this happens will depend on the coupling. We will start by considering 
an example, after which we proceed to the general solution. 


5.1 Simple examples 


Consider once again the saddle-point equation (18). To start with, assume that rn\ and m 2 


are sufficiently large so that the points -z-mi and m 2 -z he outside of the interval \-A,B\ 
(i.e. small £), where the eigenvalue density has support. This requires, in particular, that 
m 2 > 2 B and m\ > 2A In this case the saddle-point equation reduces to 


P(Z) ~-2X- 


(26) 


Normalisation then gives the condition A + B - 2A. The final condition arises from the 
integral equation, giving A - B, hence A - B - A. This solution exists for a coupling 
where the conditions m 2 > 2 B and m i > 2A are satisfied, i.e. 0 < A < ^, since m 2 < mi. 
For a larger A, the second shifted term begins to contribute in some interval, and the 
solution must therefore experience a discontinuous change. 

In general, transitions occur when one of the two shifted terms, p(-z + m 2 ) or p(-z - 
mi), are turned on. This happens when a new resonance point , reaches the interior of 
[~A,B]. Then, -z + m 2 or -z - mi coincide with -A or B. In terms of the variable -z, 
the resonance points are points located at a distance m 2 or rri\ from the endpoints of the 
interval. Physically, the phase transitions occur because beyond some critical couplings, 
extra massless particles (of masses proportional to | - z + m 2 | or to | - z - mi|) begin to 
contribute to the partition function. 

To proceed, we may consider the situation where mi > 2 A, which ensures that the first 


one of the shifted terms vanishes. The saddle-point equation (18) then reduces to 


p(U - y>(”>2 - z), 


(27) 


where we assume that p(z) is supported on the interval \-A,B], for some 0 < A < B. The 
resonance point of the leftmost interval endpoint, -A, is then given by: 


di = m 2 + A. 


(28) 
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Similarly, the resonance of the interval endpoint B is given by 


h = m 2 - B, 


( 29 ) 


but these two points will simultaneously lie inside the interior of the interval \-A,B\ 
only in one limiting case. The extremal case to have a resonance originating from -A 
corresponds to this resonance point coinciding with B, giving us a condition on m 2 as 
m 2 - B - A. When this condition is fulfilled, the resonance originating from B will be 
b l = -A. 

Thus, for m 2 < B - A, there will only be one resonance originating from the leftmost 
endpoint -A, and for 2B > m 2 > B - A, there will be one resonance originating from the 
rightmost endpoint B. [As discussed above, for m 2 > 2B , there are no resonance points 


inside the interval, and the solution is thus given by (26). 
In the case of B - A < m 2 < 2/1, we instead have: 


p(z)=- 


z e [ -A , m 2 - B ] 


a z € [m 2 - B,B] 
whereas in the case 0 < m 2 < B - A, we find: 

[ 2 A ze[m 2 + A,B\ 

Normalisation together with the integral equations allows us to fix both interval end¬ 
points in terms of m 2 ,\ and the complete expression for the eigenvalue density in the case 
where nn\ >2A is then given by: 


p ( z ) = J & z€ t~ A ’ 2m 2 - 3A] 
3 ^ z e [2m 2 - 3A, 3A - m 2 ] 


and 


p(z) = 


a z € [-m 2 , 2m 2 \ 
z e [2 m 2 , 2A] 


m 2 

— < A < m 2 


A > rn 2 


(30) 


(31) 


As seen in hgure [lj both of these cases agree well with numerics, and there is a phase 
transition at the point m 2 = A, as expected. 

Having determined A and B , we can now check the region of validity of the solution. 


For the solution (30), the condition m\ > 2 A gives the additional constraint A < mi/2. 
For the solution (31), mi > 2A requires 2() > m/3. In the complete phase diagram shown 
in fig. [3j the uniform eigenvalue density ( [26] ) is the density in the shaded region below 
the lowest blue line. The solution (30) represents the eigenvalue density in the triangular 


region above this blue line, having the green (A = mi/2) and black (A = m 2 ) lines as the 


other sides. Finally, the solution (31) is the eigenvalue density in the region above this 


black line, limited by the purple lines 2£ = m/3 and 2£ = m. 
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Figure 1: Numerical solutions for large mi with m 2 > 0. Here, ( e {0.9,0.7,0.5,0.35,0.3,0.1}™. 
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5.2 General resonance structure 


In the general case, we will have resonance points originating from both interval end¬ 
points, which may conveniently be expressed in the tables below. Table [l] contains all 
resonances originating from the leftmost interval endpoint, -A, whereas table [2] contains 
the resonances from the other endpoint, B. 


2 

-2 - 777i 

777 2 - 2 

-A 

A - 7771 

A + 777 2 

A + 777-2 

A - mi 

(-A - 2m) 

-A 

-A 

-A + 2777 

-A + 2m 

A - 2 m - 7771 

A - m i 

A - 2m - 7771 

-A + 2m 

-A + 4777 

-A + 4777 

i 

i 

A - 2m - mi 


Table 1: Table of resonance points originating from the leftmost interval endpoint, -A, 
where the first row gives the “first order” resonances, the second row the “second-order” 
etc. Normal-sized terms denote “new resonances”, whereas the smaller ones coincide with 
previous resonance points. The small one within round brackets lies outside the interval 
of support of the density. 


2 

-2 - 7771 

7772 _ Z 

B 

(~B - mi) 

7772 - B 

7772 - B 

B - 2777 

B 

B - 2 777 

m 2 - B 

-B + 2m + 777 2 

-B + 2m + 7772 

B - 4777 

B - 2m 

B - 4777 

—B + 2m + m 2 

-B + Am + 777 2 


Table 2: Table of resonance points originating from the rightmost interval endpoint, B 
(same conventions as table [I]). 

In total, the resonance points may be written as: 


\a 2 k a -i = A-m 1 -2m{k a -l) 
|a 2fca = -A + 2 k a m 

I b 2 k b -i = m 2 - B + 2 m(k b - 1) 
\hk b = B- 2 k b m 

a = 771,2 + A, 


(32) 


where k a , kb e {1,2,3,...}. From these expressions for the resonance points, it is clear 
that a and b\ cannot both lie in the interior of the interval at the same time. Rather, the 
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condition that a lies in \~A,B] is equivalent to m 2 < B - A, whereas the condition that 61 
lies inside the interval is equivalent to m 2 > B - A (i.e. exactly the condition separating 
the two non-trivial phases in the simple example considered in section 5.1), and we again 
have two separate cases to consider. We discuss these cases in the order of increasing 
difficulty; thus starting with the first one. 


5.3 m 2 < B - A 

In this case there are no resonances from the rightmost endpoint B. The only resonance 
points are then given by the a 2 k a , o, 2 k a -i, together with a. 

Here it is convenient to first use the integral equation (16) to determine B. This 
calculation is straightforward, since the argument of the sign-functions are strictly positive 
in the integration regime. Normalisation then forces all integrals to unity, and one finds: 


x<B , ->=1 by norm. 


B>-y,mi>0, -*■=! by norm. 


B = 


A rB 
2 

rri 2 <B-A , -*1 by norm. 


rti X r* 

A / dx' p{x') sign(T> - x') + — / dy p(y) sign(T> + y + m±) 
J-A 2 J-A 

= 2A, 


(33) 


A r B 

+ - 

2 


f dy p(y) sign(H + y - m 2 ) 
J-A 


in the case where m 2 < B - A. 

The interval [~A,B] will as always be divided into parts by the interior resonance 
points, which in this case are given by the points d,ci 2 k a , 02 fc a -i, where k a is limited by the 
condition that a 2 fc a -i lies within the interval. This implies k a < Let the highest 

integer which fulhls this be denoted by n, such that 


n = 


- A-C 

m 



(34) 


This means that n will be the integer number of times 2m fits in [-H,a], that is, in 
2(A - () + m. (The number of even resonances will hence be equal to n.) Define also A 
to be given by: 


A = 2(A - () + m - 2nm. 


(35) 


The total of 2 n + 1 resonances dividing the interval will be ordered as a > ci 2 n > «i > 
a 2 (n-i) > 03 > ■ ■ ■ > Ci 2 n—i > -A. Just as in the cases previously considered, the eigenvalue 
density will be piecewise constant, and we may define different constants on the different 
patches of the interval, such that: 


p{z) = 


Pk° 

Pk e 

j_ 

2A 


in the regions between an even resonance point and an odd one 
in the regions between an odd resonance point and an even one , 
between a and B 


(36) 
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Figure 2: Resonance points originating from the interval endpoint —A (shown in blue), where 
pink represents a, and the purple points are the odd resonance points whereas the magenta ones 
are the even ones. Below these are the same points under the maps z -»■ —z — rri\ and z -> m 2 - z 
respectively. The density is as such divided into a piecewise constant density, made up out of two 
sets: the density on the regions going from odd-even resonance points ( p°f7 e ), and the density 
on the regions going from the even to odd ones, denoted by p e ^°. 


as illustrated in figure [2} I 11 total, there will be n p°~ e and n + 1 p e ~ 0, s. Furthermore, the 
saddle-point equation (18) for these two different sets of constants decouple, and we are 
left with: 


2 PV + fCU + Pn + °2-k = J 
2pT + Pn-l + P°n+l-k = J 


Mk e [1, n + 1] 
Mk e [1, n], 


(37) 


where the boundary conditions are given by p^~° = p^ e = 0. 

These may be solved by a polynomial Ansatz in k, where one finds the general solution 
to be a linear function in k. The constant terms are then forced to vanish by the boundary 
conditions, giving the two sets of linear terms as: 


pV° = 


k 


pV = 


k 


(38) 


A(3 + 2 n) rk A(l + 2 n)' 

The normalisation condition on the eigenvalue density may be written in terms of 
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Pk °iP°k e i suc h that: 


71+1 71 7L 

l = A^pr° + (2m-A)E/’r-5T + l' (39) 

k =1 fc=l ZA 

Together with the definition of A, this allows us to solve for A and A in terms of m, ( 
and n, and we find: 


A = m(2n + 1) - 2((2n(n + 2) + 1) , (40) 

A = (2n + 3)(m-2C(l + 2n)), 


with B - 2X. 

Having this expression for A, we may find the critical point at which the value of n 
shifts. At the very site where n = - |, that is, we enter the region with a specific n, 

we find the following condition on 2(: 


2C 


m 

2 n + 1 


(41) 


However, it is interesting to notice that this is independent of the coupling, and so tran¬ 
sitions between different phases of this kind, with only resonance points originating from 
the leftmost interval endpoint, -A, will only occur with shifts in (/m. 

Next, consider phase transitions into a phase where the resonance points from B starts 
to move inside the interval. For fixed values of m,( (and such also n ), the condition that 
these resonances will remain on the outside of the interval corresponds to m - 2( < B - A, 
which in terms of the coupling may be written as: 


(n + 1 )(m - 2( t (n + 1)) < A. 


(42) 


Therefore, as A decreases, we leave this regime and enter the next one, where resonances 
appear from both endpoints. 

The general eigenvalue density in this case, for an Fl-parameter in C < y an d coupling 
satisfying equation (42), will be given by: 


P(z) 


A(3+2 n) 
k 

A(l+2 n) 

J_ 

„ 2A 


VZ 6 [a 2 (fc-l), 02(n-fc)+l] 

V z e [a 2 ( n -A:)+i,02fc] 

Mz e [a,2A] 


k e [1 , n + 1] 

ke[l , n\ , (43) 


where ao and a_i 
equation (34). 


should be interpreted as -A and a respectively, and n is defined by 


5.4 777-2 > B - A 

This situation is more complicated. Although a now lies outside the interval, there are 
now interior resonance points originating from the rightmost interval endpoint B as well. 
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The complete set of resonances in this case is given by the expressions in equation (32), 
namely: 


\a 2 k a - 1 = A - mi - 2m(k a - 1) V 0 < fc a < + \ 

I a 2ka = -A + 2 k a m VO <k a < 


jb 2 k b -i = m 2 - B + 2 m(k b - 1) 
\b 2 k b = B- 2 k b m 


V 0<fc 6 <^ + i 


(44) 


where k a ,k b e 1,2 ,... and the upper limits comes from requiring the resonances to lie 
inside the interval of support of the eigenvalue density \-A,B\. 

When m 2 > B - A, it is easy to see that 


A+B A-C 1 B + C 1 

- < - + - < - + -, 

2 m m2 m2 


(45) 


and both k a ,kb are then limited by 4^. We again define n as the integer part of this 
number, 


n = 


'A + B' 
2m 


(46) 


such that a 2n ,b 2n denotes the final even resonances. However, in some cases, the following 
odd resonances may lie within the interval [~A,B] as well. This gives rise to three different 
cases that are examined in appendix [A] in detail. These are: 


. Case I (n = ], 


WC i ; b+c 

m 2 m 


< n ): 


A =A 
B =A 


2n + 1 


n +1 
2 n + 1 


n + 1 


-2C n 

+ 2 (n 


(47) 


Case II (n = , 


a-c 

m 


\<n< — 




. . 2n + 1 . 

A -\ -2 Cn 

n + 1 

B - - m + 2((n + 1) + A 


2n + 3 
n + 1 


(48) 


Case III ( 


n 


= \a±b] 

L 2 m i 


n < 


A-C 1 ' B+Q lx. 
m 2 m 2 /' 


A = - m - 2((n + 1) + A 
B = - m + 2£(n + 1) + A 


2n + 3 


n + 1 
2 n + 3 


n + 1 


( 49 ) 
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5.5 Phase transitions 


The condition m 2 > B - A also manifests differently in the three different cases. In cases I 
and III, this turns out to be completely independent of the coupling, whereas this is not 
the situation in case II. However, this is to be expected since the second situation is quite 
similar to the situation with resonances only from the interval endpoint -A, considered 


in section 5.3 The condition found upon the coupling in case II is precisely 


A< \t 


-V 4 = (n + l)(m-2C(re+ 1)) , 


(50) 


which is precisely opposite to the condition (42) in section [573 Thus, as the coupling 
decreases for a fixed m,( and, therefore, n, eventually, the inequality is saturated in (50), 


the resonances from B move inside the interval and a phase transition occurs, leading to 
case II of this section. 

The conditions for the two other cases considered in this section are 


Case I: 


Case III: 


c< 
c < 


m 


2{n + 1) 
m 

2{2n + 3) 


(51) 


In this situation, where ( < y, and m 2 > B-A, there are clearly phase transitions. This 
is expected, since they do appear in the situation with vanishing Fl-parameter. However, 
there are different kinds of phase transitions: both in between the three different cases 


(described in detail in sections A.l - A.3), within the same value of n, and one where the 
value of n changes. 

Let us start with some fixed value of n, such as n - [4^y], and let us consider the 
case where n is larger than - |. Then, as the coupling increases, so does | 

(growing linearly with coupling) for some fixed n. At some point this will surpass this n, 
and a phase transition occurs that will take us to case II above. This will happen for the 
coupling: 


>Lir j( n+ 1 )(™- 2 C), 


(52) 


which in the case of vanishing Fl-parameter simply corresponds to the situation where 
another factor of m fits inside the interval. 

From there onwards, as the coupling grows further, so will the factor | (also 

growing linearly with A), and at the point 


A c 


n-.m= 2 ^ n + 1 )( m + 2 f), 


(53) 


another phase transition occurs, taking us into case III. It is worth noticing that this phase 
transition does not occur in the case of vanishing Fl-parameter, but rather coincides with 
the one between case I and II, simply because the second case, where n £ -1], 

never occurs for ( = 0, since this interval then is empty. 
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The final phase transition then occurs as one goes from case III back to case I, 4^ 
has grown to the point that the integer part of it changes, and that n increases with one. 
This happens at the point = n + 1, occurring at 


A 


C 

III-+I 


m(n + l)(n + 2) 
2 n + 3 


(54) 


This corresponds to the result obtained in (T| for the transition from m fitting a total 
of 2n + 1 times inside the interval [-A,A] to m fitting a total of 2n + 2 times, for £ = 0. 
(Again, in this special case, this corresponds to another multiple of m fitting inside the 
interval.) 

In order to determine the order of these phase transitions, we need to study the analytic 
properties of the free energy, 


F 


1 

RN^N ~ 2 


In Z, 


(55) 


at the critical values of the coupling. It turns out that both the first- and second-order 
derivative of F are continuous, whereas the third-order derivative is not. Let us denote 
by A F the difference between the subcritical and supercritical free energy in each phase 
transition. For the phase transitions occurring between the case with only resonances from 
- A , for some given n, (presented in section 5.3) to case II with this same n presented 
above, occurring at A^ = (n + 1 )(m -2 ((n+ 1)), we find: 


dx AF 


aA 


= d'l AF 


aA 


- 0 , 


(56) 


f)3 A Tp _ _ 

A AA (n+l) 4 (m-2C(n + l)) 2 ' 

Note that the apparent singularity at m - 2 ((n + 1) is outside the region where these 
solutions apply, m/(2n + 3) < 2£ < m/(2n + 1). 

Similarly, at the critical points between the cases I, II and III above (described by 


equations (52) - (54)), we once again find a discontinuity at the third derivative, whereas 


all lower derivatives are continuous: 


dl A F 


32 


a aJ(")-h(") (n + l) 3 (m - 2() 2 

32 


(57) 


dl A F 


A“ ( n + l)3( m + 2 £)2 

2(2n + 3) 5 


dl A F 


A in(n)-»i(n+i) m 2 {n + l) 4 (n + 2) 4 


Therefore, just as in the case with vanishing Fl-parameter, these phase transitions are all 
of third order. 
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Finally, one may also look at transitions between solutions (43) of section 5.3 with 


different n, as the FI parameter ( is increased along lines of constant A. The solution (31) 


represents the case n - 0 and is valid in the region 2£ > m/3, A > m 2 . The case n - 1 can 
be readily found from the general formulas of section 5.3[ We obtain 


p(z) = 


5A 

J_ 

3A 

_ 2 _ 

5A 

J_ 

2A 


Z € [ — A 1 fll ] 
z e [ai, a 2 \ 
z e [a 2 , m 2 + A] 
z € [m 2 + A, 2A] 


A > 2m - 8C , — < 2C < — , 
5 s 3 


(58) 


with 


A = 3m - 14£ , ai = 2m - 16C , a 2 = 14£ - m 


(59) 


In crossing the line 2( = m/3 from lower to higher values of (, the solution changes from 
(58) to (31) (in fig. [3j this corresponds to crossing the purple line at 2( = m/3 that 


begins at an hexa-critical point). The discontinuity in the free energy shows up in the 
third derivative with respect to (. Again, for convenience we first compute d\F which 
has a simple local expression, given by ~(z 2 )/ A 2 . We find that first and second order ( 
derivatives are continuous, whereas 


d 3 c (A d x F) 


2C=? 


432 

W" 


(60) 


More generally, for the transition from the n - 1 ->• n solutions (43) we find 

16(2n +l) 3 


df(Ad x F) 


A 3 


(61) 


2n+l 

Thus the quantum phase transitions between these phases are also of the third order. 


6 Summary of the Analytically Continued Model 

6.1 Case Ai = A 2 

We have herein seen that the phase structure found in [l] is significantly enriched when 
the theory is deformed by a Fl-parameter, (. The saddle-point equations then loose the 
reflection symmetry present at £ = 0, and the behavior in the decompactihcation limit is 
highly dependent on the value of this new parameter. As £ 0, we recover the results 

of 0. 

For an Fl-parameter large enough, ( > the theory is free from phase transitions. 
However, for |£| < phase transitions appear. The precise appearance of the eigenvalue 
density and the position of the phase transitions in phase space depend on A and ( in 
relation to m. 
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247m 


Figure 3: Phase diagram of the analytically continued model. Black critical lines separate 
phases of case II (below) from phases with only resonances from -A (above). The shaded 
region corresponds to the constant eigenvalue density (26). Crossing the blue, green and red 
lines corresponds to phase transitions from cases with resonances from both endpoints, I -*■ II, 
II -*■ III and III -*■ I (where n increases by one), described in section 5.4 The purple dots are 


hexa-critical points. The purple (vertical) lines separate phases described in section |53| differing 
in one unit in the value of n, representing solutions with 2n + 1 resonances from the leftmost 
interval endpoint. At the orange triangles the eigenvalue density has been computed numerically 
and compared to theory (see fig. [5j. 


To illustrate this dependence, one may consider a phase diagram with the dimension¬ 
less axes A/m and 2£/m (figure [3]). Phase transitions then occur on certain critical lines 
in this phase space. There are no phase transitions in the region where the Fl-parameter 
satisfies ( > y, and so there is no use to show the phase diagram further than 2C/m = 1. 


For 2£/m < 1, the only phase with no resonances is the constant eigenvalue density (26), 
occurring in the shaded region of fig. [3j At the purple dots, there are six coexisting 
phases: they are hexa-critical points. 

It is worth noting that the phase transitions between cases I, II and III only occur in 
a specific order (see fig. [4]). The number of phase transitions undergoing as the coupling 
A runs from 0 to infinity depends on the value of 2C/m, and the maximal n max that will 
occur is given by the largest n fulfilling < A^(n), which, by using equations (|50| 


and (54), implies: 


n 


m 


AC 2 . 


(62) 


The number of phases is then given by 3(n max + 1), It approaches infinity as 2 C/m -*■ 0. 
This is to be expected in order to match the £ = 0 case [l]. However, for any non- 
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Case / 

Case III <- Case II 

V_ Only 
resonances 
from -A 



Figure 4: Sequences of phase transitions in the theory at constant (. As A grows, the cycle 
continues until the value of n max where, for given m,(, one has A* I ~ > " 4 (n max ) < Aj I_>III (n max ), 
after which the cycle ends. 


vanishing Fl-parameter, there is only a finite number of phase transitions in the theory 
(and no phase transition for 2( > m). 

The grey, vertical line in figure [3] at 2 (/m = 0.22 illustrates an example on the phase 
structure of the model, and how to read the diagram. For 2 (/m = 0.22, one has n max = 1 by 
(62), hence 6 different phases along the grey line. Starting at the bottom of the diagram, 
we are in the trivial phase (case I with n - 0). As the coupling grows we moves upwards 
along the grey line until we cross the first blue line. This corresponds to a phase transition 
into case II with n - 0. After crossing the green line, we enter III with still n - 0, and 
when crossing the red line, we move back into case I, but now with n = 1. Crossing the 
next blue line takes us to case II. n= 1. However, instead of crossing another green (and 
thereafter red) line, we rather cross a black line. Beyond this point, all resonances from the 
rightmost interval endpoint move outside the interval of support of the eigenvalue density, 
and we thus get into the phase described in section [A3] with n = 1, having resonances only 
from the leftmost interval endpoint -A. 

As the coupling grows further, no new resonance points enter the interval, and the 
system never leaves this phase. The orange horisontal lines in figure [3j together with the 
orange triangles, represents points at which numerical calculations have been made in 
order to compare with theoretical calculations. The eigenvalue density for the six phases 
present for values of m, ( such that 2 (/m = 0.22 are presented in figure [5j and we find 
an excellent agreement with our theoretically derived densities in sections |5.3 5.4 and 
appendix [A] The smoothness of the curves visible in figure [5] has to do with finite-size 
effects. 

As a check, using equations (47)-(49), we may compute the values of n in these different 
phases, and these results are summarised in Table [3] The results follow the expected 
pattern, where n changes in going from case III to case I. 

The results of section 5.4 (and appendix [A]) explain some of the peculiar properties of 
this model in the special case of ( - 0 : as resonances from both the left- and right interval 
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A/m 

0.22 

0.5 

0.63 

0.73 

0.96 

1.2 

n 

[0.22] = 0 

[0.61] = 0 

[0.89] = 0 

[1.095] = 1 

[1.53] - 1 

[1.85] = 1 

Phase 

ln=0 

Hn=0 

IHn=0 

In=i 

IIn=l 

(only A-resonances) n= i 


Tabic 3: Sequence of transitions and corresponding n = [4^] in cases I - III and n = 
[tT ~ l] the phase with only resonances from one endpoint. 

endpoints appear in the interior of the interval, the interval is divided into four quali¬ 
tatively different regimes: two in between resonances originating from the same interval 
endpoint, and two in between resonances originating from different interval endpoints. 
As the Fl-parameter goes to zero, these regimes become pairwise indistinguishable from 
one another due to the presence of reflection symmetry. The different origins of these 
regimes are indeed visible in the eigenvalue density, as “odd” and “even” patches of the 
eigenvalue density behave significantly different from one another; in addition, the criti¬ 
cal couplings at which phase transitions occur are different for “odd” and “even” phases 
(where mass parameter m fits inside the interval an even- or an odd number of times). 
This is all explained here by the fact that these regimes actually originate from different 
sets resonances; some originating from the leftmost endpoint, and some from the right¬ 
most endpoint of the interval. In this way, the general solution with ( + 0 explains a 
phenomenon which appears to have no deep reason in the case of vanishing Fl-parameter. 

Furthermore, one may note that the green and blue lines in figure [3] coincide when 
( = 0. Here, the only phase transitions occurring are the ones between the cases I and 
III. This is required to agree with previous results for ( - 0 [l]. 
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Figure 5: Eigenvalue densities determined numerically for the six phases present with 2£/m = 
0.22 and A/m e {0.22,0.5,0.62,0.73,0.96,1.2}, compared with eigenvalue density determined 
analytically, for the cases I n =o,IIn=o,III n= o,In=i,IIn=i)^4n=i (grey lines). [The smooth edges of 
the eigenvalue density are due to the finite radius R used in the calculations. Here, mR = 50. 
For higher mR the numeric is unstable.] 


6.2 Case Ai ± 


6.2.1 First analytic continuation 


It is interesting to ask what happens in the case when the couplings are different. We 
assume Ai, A 2 are real (as discussed, for general complex values of Ai, A 2 the saddle-point 
equations become complicated to solve, even numerically). Then we have to compute the 
partition function in the convergence region where Ai. 2 > 0. We further assume, with no 
loss of generality, Ai > A 2 . 

First, let us consider the limit of very large Fl-parameter. In the case of ( -*■ oo, both 
shifted terms of the saddle-point equations 0 vanish, and we are left simply with 


Px — 




(63) 


Using the integral equations (16) with x - B,y 


D, together with the normalisation 
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conditions for the densities, one finds the interval endpoints to be given by: 


B - D - \i + X 2 (64) 

A = -C = Ai-A 2 


Consider what happens when the Fl-parameter decreases, 
where the above solution (63) may cease to be valid: 


Then there are two points 


• A resonance point moves inside the region of support for the eigenvalue density. 
The first point to do this is the one originating from the leftmost interval endpoint 
of the other interval. That is, A-m 2 (or C -m 2 ) moves inside the interval [-C,0] 
(or equivalently [-A,T>]). 


The sign-functions change inside the integration regime in (16). 


Given the solution in the case ( -> 00 , we find that the first one of these situations 
occur precisely when m 2 < 0 and, just as in the case with equal couplings, there will be a 
phase transition in the theory as soon as ( < y. 

As for the second possible point of failure for the solution (63) for large (, one may 
easily show that the sign-functions will not change signs inside the integration regimes 
until C < y - A 2 , C ^ y-Ai, which for positive couplings always happen after the phase 
transition originating from resonance points entering the intervals. This shows that the 
theory with ( > y does not have phase transitions, generalizing the result found in section 
[4] to the case of arbitrary Ai, A 2 > 0. 

We now derive the solution in this new phase, where ( decreases just below y. Assume 
C = y - u, for some u, satisfying 0 < u < y (together with some other conditions which 
will shortly be specified). Then one finds 


mi =2 m - 2 u (> 0) (65) 

m 2 =2 u. 


(It is here worth noting that both terms shifted by mi in the saddle-point equations (16) 
will vanish, at least for small u, since we know A = -C in the limiting case of u = 0.) 

We wish to consider a situation where the resonance point from the leftmost interval 
endpoint of one interval lies in the interior of the other one (i.e. A + m 2 < D , C + m 2 < B), 
which gives us the additional conditions on u mentioned above. This is equivalent to 
saying that the resonance points originating from the rightmost endpoints of the intervals 
will lie outside the region of support for the densities. As long as these conditions are 
satisfied, the integral equations may once more be used straight away to obtain 


B — D — Ai + X 2 


( 66 ) 


just as in the case of ( > y. 
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The intervals \-A,B],[-C,D] will then be divided into two parts each, and using a 
piecewise constant Ansatz for the eigenvalue density on these patches, together with the 
normalisation condition, one obtains: 

M z e [—Ai + A 2 — 2 u , — Ai + A 2 + 4m] 

Vz e [-Ai + A 2 + 4m , Ai + A 2 ] 

Vz e [—A 2 + Ai — 2 u , — A 2 + Ai + 4 m] 
yz € [—A2 + Ax + 4 u , Ax + A2] 

Thus, as the Fl-parameter decreases just below y, a phase transition occurs, even 
though Ai ± A 2 . Additional phase transitions are undergone as resonance points continue 
moving inside- or outside of the region of support of the eigenvalue densities. The quali¬ 
tative behaviour of the theory is therefore similar to the Ai = A 2 case studied in previous 
sections. 



6.2.2 Second analytic continuation 


Consider for simplicity mi = m 2 = m, i.e. £ = 0. The saddle-point equations (19) imply 
the following functional equations: 

: + m) + -/ 

x 1 

z + m) + -/ 

The trivial solution is 

1 


pM= 2 y + ? 


Px(x) = 


X 6 (-Ox, Ox) 


2ax ’ 

Py(y) = o — ) 1/6 (-a 2 ,a 2 ) 

Z«2 

and holds provided ai + ct 2 < m. 

When ccx + «2 > m, the solution has two patches: 

(A + A x€[a,A] 


m) , 

(67) 

m). 

(68) 


(69) 


IA + A 


(70) 


a-m-B 
b-m - A 


cki (2m - a 2 ) 

2aq + e»2 

ct2 (2m - cti) 

Ox 2«2 


«2 (aq + m) 
2 ax 6 X 2 

CXx («2 + TTl) 

ot\ + 2cx 2 


( 71 ) 

(72) 
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where we only exhibited the region x > 0, since in this ( - 0 case the density has reflection 
symmetry. 

Thus there is a phase transition for generic an, 012 > 0. Note that the solution requires 
ai, a 2 to be positive, as the eigenvalue densities must be positive. The free energy can 
be computed from the formulas 

d « F -^' < 73 > 


One can check that the resulting expressions satisfy the integrability condition d ai d a2 F = 
d a2 d ai F. For the phase with uniform density (69), we find 


Fi - —(a 1 + 02 ) 


(74) 


For the phase (70) 


Fn = 


-6 (ai + 02 ) m? - 3aia2 m + aia2 (ai + 02 ) + 2 m 3 
12 ( 2 ai + 012') (a 1 + 2a2) 


(75) 


It follows that the first and second derivative are continuous at the transition point 
ai + a 2 = m , whereas the third derivative has a jump: 


OlAF 


ct2=m-ai (2/n - ai)(m + ai) 


AF = Fi~ Fa 


(76) 


Therefore the transition is, as in the first analytic continuation, third order. It occurs for 
generic, positive values of ai, a 2 at ai + a 2 = m. Note that this phase transition is not 
meaningful in the ABJM case: it occurs when ai + 012 > m , which is never the case in 
ABJM where ai + a 2 = 0. 
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7 Mass/FI-deformed ABJM theory does not have 
phase transitions 

In previous sections we solved the mass-deformed ABJM model with a non-vanishing FI 
term at large N , with the couplings analytically continued into the complex plane. This 
is a standard approach in studying the large N behaviour in Chern-Simons theories (see, 
e.g., (4|[l8[25j,26|). We have found a rich structure of phase transitions, but an important 
question is whether the ABJM model with physical couplings and masses exhibits phase 
transitions as the coupling is varied. For the mass-deformed models, analytic continuation 
back to physical couplings is not straightforward, due to the existence of poles originating 
from the hyperbolic functions in the partition function. Clearly, it would be more desirable 
to have a direct solution of the large N ABJM model with the original parameters, with 
k\ = -k 2 integers. Solving the saddle-point equations in this case turns out to be very 
complicated, because eigenvalues seem to be distributed in cuts in the complex plane, 
with non-uniform N dependence. In this section we will argue that physical ABJM 
theory deformed by arbitrary mass and FI terms is free from phase transitions. 


7.1 General case 


We start with (J2]), with m < 2(. In section 6.2.1 we have shown that this theory does 
not have phase transitions for generic Ai, A 2 , real and positive. Therefore the free energy 
is a smooth function of the couplings Ai,A 2 , indicating that there should not be phase 
transitions after Ai j2 -*■ e ±zlp Xi )2 analytic continuation back to the physical ABJM model 
when m < 2£. But the partition function of the ABJM model Zabjm( 2C i m ) (with N\ - N 2 ) 


has the symmetry (10). The symmetry holds for any N , in particular, in the planar limit 
N -»• 00 with fixed k/N. If for given ( and m, Zj^ B j M (2(,m) is a smooth function of 
k/N, then so is Z^jM^m^O. This indicates that there should be no phase transitions 
in ABJM neither in the opposite regime when m > 2(. 


7.2 Massive U( 2) x U( 2) ABJM 

We now consider a case where the partition function can be computed exactly, namely the 
N - 2 case, i.e. U (2) x U (2) ABJM deformed by mass and FI terms. We will first consider 
the analytically continued model with real, equal couplings, particularised at N = 2. This 
corresponds to the second analytic continuation in the Chern-Simons levels, defined in 


section 3.2 and further studied in section 6.2.2 (now the first analytic continuation cannot 
be used because the ranks of the two gauge groups are fixed from the beginning). We will 
show that, in the same decompactihcation limit, the model also exhibits phase transitions 
of the same nature as the large N model. Then we will discuss the U (2) x U (2) ABJM with 
physical coupling and parameter deformations, and show that there is no phase transition 
in this case. 
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7.2.1 The analytically continued model with equal, real couplings 

Consider the gauge group U( 2) x U( 2), and analytic continuation in the Chern-Simon 
levels to equal, real couplings. The partition function ([Tj) takes the form 


z=i 


d 2 ii d 2 v sinh 2 f (/xi - fx 2 ) sinh 2 f (iq - z/ 2 ) e ^ ) 

4 J (27t) 2 (27t) 2 riij=i cosh(f (fj,i - v 3 + m)) cosh(f (/^ - - m)) 


/ 


(77) 


where the /^-dependence has been restored. In this subsection, for simplicity we set the 
Fl-parameter to zero, since this case already illustrates the main point. This is the analog 
of the models with cq = a 2 studied at large N in section 4 of (T]), and briefly reviewed in 
section 3.2 They have similar phase transitions as the Ai = A 2 models obtained by the 
analytic continuation in the gauge group ranks, discussed in detail in sections 4,5, 6. 

Although phase transitions in matrix models typically arise at large N, this type of 
phase transitions due to the contribution of extra massless multiplets at certain couplings 
have also shown up in some finite N examples 11 |l4|. The reason is that the contribution 
of an extra massless multiplet produces a singular behaviour at any N, even when N - 2. 
We now wish to see if the present model also has phase transitions. 

The non-analytic behaviour arises upon taking a suitable decompactihcation limit R -*■ 
oo, where the integral defining the partition function is dominated by large expectation 
values /ij, zq. In this limit, the hyperbolic cosine functions in the denominator get replaced 
by the non-analytic functions |exp|/q - Vj ± m\, which produce non-analytic behaviour 
when the coupling is such that /i,; - Vj at the saddle-point hit ±m. We therefore assume 
the same scaling we used in the large N case: 


g = XR , 

R -> oo , A fixed. 

For large R , Z may be written as: 

7-4 f d 2u -Rjj/Jt-H+ht-H-! sL=i EA/h^)] 

J (27t) 2 (27t) 2 

and the integral is dominated by a saddle-point. The saddle-point equations are 


(78) 


A 


1 1 2 

/q = sign(/q ~H2 ) - - ^ [sign(/q - v 3 -m)+ sign(/q - v 3 + m)\ 

Z 3 =1 

1 1 2 

t h2 = -sign(/ii -/i 2 ) - - X! [ si S n (h2 - Vj ~ m) + sign(/i 2 - v 3 + m)] 
A 2 j=1 

1 1 2 

Th = sign(z/i - v 2 ) - - Y, [sign(z/i - \i 3 - m) + sign(/y, - /q + m)\ 

A 2 j=1 

1 1 2 

y v-2 = -sign(i/i - v 2 ) - - Y, [sign(z / 2 - [X 3 -m) + sign(z / 2 - ^ + m)\ 
A 2 j=1 


(79) 

(80) 
(81) 
(82) 
(83) 
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We can assume with no loss of generality Hi > /i 2 and V\ > u 2 . If the eigenvalues are 
sufficiently small, their difference will be less than m. Then the sign functions containing 
m in the argument cancel out and we find the solution 


Aii = A, H 2 = -A , vi = A , u 2 = —A . (84) 

Thus this solution holds foi0 

A<y- (85) 

When A > y, the solution will change, because the difference of eigenvalues can overcome 
m and in this case the sign functions will contribute. In this A > y regime, the absolute 
minimum of the action is given by 

771 711 

= = V2 = v 2 = -—. (86) 

In this case the arguments of some sign functions vanish. The action is not differentiable 
at this point and the minimum must be found by inspection. 

Therefore the theory contains two phases. The free energy F = -j, hi Z in each phase 
is given by the action evaluated at the minimum of the potential. We find 

_ |4m - 2A A<f 
\2m + ”£ A>f 

This implies a discontinuity in the second derivative 


A F 


= d\AF 


\ _ m 
2 


\ m 
A “ 2 


-0 , 


d 2 x A F 


\_m 

2 


m 


(87) 


Therefore, we conclude that the analytically continued 1/(2) x[/(2) mass-deformed ABJM 
model presents phase transitions. 


7.2.2 U( 2) x U( 2) ABJM model with physical couplings 

The partition function of the U( 2) x U( 2) ABJM deformed by mass and a FI term com¬ 
puted by localisation is given (see ([l]) with N - 2) 




/ 


d 2 /i dV 


sinh 2 sinh 2 ^ 


\2 nr2 ( fH-Uj+m J cqg1i (' Vi-Vj-m j 


(EiW+Eii'-) 


g 47T ■ 


4J (2tt) 2 (27t) 2 riij=i cosh(- 2 , 2 

Using (|8]) for A = 2, the partition function can be written in the following form 

zF( Zl -z 2 ) , 


( 88 ) 


(89) 


*For two independent couplings Ai, A 2 , the solution is Hi = -/r 2 = Ai, V\ = -v 2 = A 2 , with Ai + A 2 < to, 
which is the N = 2 analog of the solution (691. 
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with 


Zi = J dr\ dT‘> 


g-ifcmj(ri+T2) 


and 

^2 = J dr | (It 2 


cosli(7rfcri) cosh(7r kr 2 ) cosh 2 (^) 

g-ifcm2(Ti+T2) 


(90) 


cosh(7rfcr[) cosh(7rfcT 2 ) cosh (tt(ti - r 2 ) - v 2 -)') cosh (^(t! - r 2 ) + 

(91) 


mi = m + 2£ , m 2 = m- 2( 


(92) 


It is important to note that the derivation that leads to this form of the partition function 
holds if and only if the Chern-Simons levels of both gauge groups in U( 2)^ x U(2)k 2 are 
opposite, i.e. k 2 = -k\. 

The first integral can be computed by using the Fourier integral (|6|). We get 


Zi = 


k 2 cosh 2 (^) cosh 2 (^f) ' 


(93) 


Next, consider the calculation of Z 2 . Note that we have reduced the original four integrals 


to only two integrals T\ and r 2 . The integral (91) can be carried out by defining variables 
u = T\ + t 2 , v = t i - r 2 . Then the integral over u is a Fourier transform than can be 
computed explicitly: 


/ 


du 


g-ikrri2U 


2 sin (km 2 v) 


cosh 1 y(u + v) cosh ^(u - v) ksmh^kv) sinhm 2 


Hence 


Z 2 — 


/ 


dv- 


sin (km 2 v) 


k sinhm 2 J sinh(7r kv) cosh(7rn - cosh( 7 ™ + 

Restoring the R dependence, and rescaling v -*• Rv, we have 
2A K 2 r , sin(m 2 uR/A) 


Zo = 


/ 


dv 


sinh (m 2 R) J sinh(7rn/A) cosh(R(nv - ^)) cosh(R(nv + 


(94) 


(95) 


(96) 


where we defined A = l/(kR) = fixed. At large A, this integral is not dominated by 
a saddle-point; nevertheless, it can be computed exactly. For R -»• 00 , the product of 
the hyperbolic cosines in the denominator becomes proportional to a step-function with 
support in the interval (-5=)y,?y). The sin(7rn/A) in the denominator can be replaced by 
ttv /A (as can be seen e.g. by a change of integration variable v -*■ x/R). The resulting 
integral can be carried out explicitly, with the result 


qo h 

Z 2 = e --R(KW"»al) Si( — |mim 2 |R 2 ) , 

7r k 2 2n 


(97) 
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where Si(z) is the Sine integral function. Combining with Z\ at large R , we finally obtain 
Z = A e -tf(l"*i|+|"w|) 1 1 - —Si( —|m 1 m 2 |i? 2 )') 

fc 2 \ 7T 27T / 

QO b 

» e _R(|mi|+|m2|) cos( —|-mim 2 |i? 2 ) , (98) 

/C|mim 2 |/t 2 Z7r 

where we have used the asymptotic expansion of the Sine integral function. 

Thus we have evaluated the U(2)k x C/(2)_^ ABJM partition function with both mass 
and Fl-parameter deformations in the large /? limit. In particular, this shows that F = 
- In Z is an analytic function of the coupling k (or A). Therefore, for physical couplings, 
the model does not exhibit phase transitions. More generally, since the Sine integral 
function is an entire function in the whole complex plane, the theory does not have phase 
transition in any region of the /e-complex plane. Thus, phase transitions seem to be 
absent in U( 2) x U(2) ABJM theory, where the Chern-Simons levels of the gauge group 
are opposite integers. 


Remark: The partition function (98) has zeroes, which, at large R, are located at 


kmim 2 R 2 = 


?7iim 2 i? 

A 


7r 2 (2n + 1) , neZ 


(99) 


These are Lee-Yang singularities])] It is easy to see that the partition function has zeroes 
also at finite R. It would be interesting to get further insights on their physical meaning. 
They appear to be resonances occurring at special values of mR , perhaps associated with 
Kaluza-Klcin excitations in § 3 . They seem to arise by virtue of the fact that: a) Chern- 
Simons theory has imaginary coupling g = 2iTi/k] b) the theory has mass/FI deformation; 
c) the theory is on a compact space. As a small check, one can see that the partition 
function of other Chern-Simons theories with massive matter also exhibit similar zeroes. 


In particular, for U{2) Chern-Simons theory with fundamental matter one finds 14 


7 u( 2 ) 

Zj {k= 2 ) 


„ r, r, p / im 2 R 2 \ ( im 2 R 2 A 

87i 2 e 2mR I e 2 tt - 1JI e ^ + 1 1 

(e 2 mR - l) 2 


( 100 ) 


which, indeed, has an infinite number of zeroes at m 2 R 2 = 7r 2 n, m 2 R 2 = 7r 2 (4n - 1), 
n = 1 , 2 ,... 


8 Concluding remarks 

In the first part of this paper we have studied the general solution to functional equations 
of the form 

p{z) - \p{~z ~ mi) - \p(~ z + ™ 2 ). (101) 

^We thank K. Zarembo for this remark. 
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where p(z) is a unit-normalised density supported on some interval \-A,B] along the real 
axis. We have shown that this equation describes the large N limit of the mass- and 
FI-deformed ABJ partition function Z( Ai,A 2 ) with couplings Ai = 2niNi/k, X 2 = 2niN‘2/k 
analytically continued to the complex plane, in the particular region where Ai = A 2 e R. 
The study herein generalises the discussion of [I], corresponding to the case mi = m 2 in 
this notation, and this generalisation turns out to be highly non-trivial since it amounts 
to giving up the reflection symmetry around the origin which drastically simplified the 
analysis in the special case of m\ = m 2 . 

Unlike the case of vanishing Fl-parameter, for ( ± 0 the theory exhibits a finite number 
of phases as the coupling is increased from 0 to infinity. The structure of these phases is 
summarised in section [6] and in the complete phase diagram given in figure [3] One result 
obtained here of particular interest is that the FI deformation alone, while it introduces a 
mass scale, does not generate phase transitions. Furthermore, in section [6j we have also 
considered the cases of generic (real and positive) couplings Ai, A 2 , as well as the case of 
analytic continuation in the Chern-Simons levels instead, and showed that the qualitative 
picture is similar to the equal coupling case of previous sections. 

An important question concerns the implications of these results for the physical 
ABJM theory, where the couplings are opposite (Ai = -A 2 ) and purely imaginary. We 
have presented an argument showing that there are no phase transitions in this case. The 
argument relies on the symmetry of the original ABJM partition function under exchange 
of FI- and mass-deformations, along with the fact, derived in section 6.2.1 , that the ana¬ 
lytically continued theory is free from phase transitions when ( > m/2. To further clarify 
this issue, we have computed the partition function for the gauge group U( 2) x U( 2) ex¬ 
actly. While the theory has phase transitions of the same nature as the large N model in 
a region of parameter space (in particular, at equal real, couplings), there are nevertheless 
no phase transitions in the physical ABJM case, where the Chern-Simons levels are k and 
-k. 

The U(2)k x U( 2)_fc example explains why the U(N)j- x ABJM model should 

be free from phase transitions. It would be interesting to get further insights on the whole 
picture and on the analytic structure of the mass/FI-deformed ABJM partition function. 
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A The phases for m 2 > B - A 


In section [5~4| we found that, when m 2 > B - A, even and odd resonances can get into the 
interval \-A,B] of support of the density in three different ways: 


• Case I: Both points a 2n +i and b 2n+ i lie outside the interval, which is equivalent to 

( 102 ) 


A - ( 1 , B + ( 1 ^ 

-<-< n. 

m2 m2 


Case II: The point a 2n+ i lies outside the interval, but b 2n+ i lies inside the interval, 
corresponding to 

A -Cl „ B + C 1 

2 m2 


m 


(103) 


Case III: Both a 2n+ i and b 2n +i lie inside the interval, that is, 

1 

2 ' 


, A~C 1 , B + C 

n <-< - 

m2 m 


(104) 


In what follows we examine each case separately. 


A.l Case I: n = [^], 


AC _ I. < Ah _ I < n . 

m 2 m 2 


In this case, the ordering of the resonance points in the interior of the interval will be 
given by: 

-A < b 2n < b\ < a 2n -i < (105) 

< a 2 < & 2 (n-l) < b 3 < a 2 ( n -i)-i < • • • < 

< « 2 n -2 < b 2 < b 2n -1 < a\ < a 2n < B, 

which is illustrated in figure [6] With a piecewise constant Ansatz for the density on the 
form: 


p(z) = 


a/3k 

Vz e [a 2 fc,fe 2 (n-fc)] 

k e [0, n] 

2k 

€ \b 2 ( n+ i-k),b 2 k-i\ 

k e [l,n] 

2otk 

Vz e [52fc-l,02(n-fc)+l] 

k € [l,n] 

<*k 

\/ Z e [a2(n-fc)+l) a 2fe]) 

fc e [l,n] 


(106) 


with -A - a 0 ,B - b 0 , the saddle-point equation (18) for a/3k,2k,2 a k,ctk takes the form: 

1 


2>oif3k + (3oi n —k + /! cy 


n+l-fc — 


A 

2 2k + 2n+l-k + 2n+2-k = 

2/3cifc + a(3 n -k + cx2n+i-k — — 

2«fc ^n-k ®n+1 -k ~ ^ 


k e [0,n] 
k e [l,n] 
k e [l,n] 
k € [l,n], 


(107) 
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Figure 6: Resonance points originating from the interval endpoint -A and B (which are shown 
in blue). Resonances originating from —A are shown in magenta (even, a 2k ), and purple (odd 
resonances, a 2 k-i)- Similarly, resonances originating from the rightmost interval endpoint, B, 
are shown in yellow (even, b 2k ) and red (odd, & 2 fc-i)- These resonances exist for all k < 
Furthermore, the resonances under the maps z -> -z-m i and 2 -+ m 2 -z are shown underneath. 


with the boundary conditions /3 n+ i = a 0 = (3a n+ i = /3a 0 = 0. The decoupled equations, 
for a, (3, are straightforward to solve, and from the two remaining coupled equations, one 
finds that the general solution for f3a k is a linear function in k. Enforcing the boundary 
conditions, both integration constants vanish, completely specifying (3a k , and thus also 
a/3k- Altogether, one then Ends the eigenvalue density as: 



1 

2A 

Vz e [a 2 k,b 2 (n-k)] 

k e [0,n] 

Pk 

_ 1+n-k 

A(l+2n) 

Vz € [&2(n+l-fe) ,^2fc-l ] 

k e [1 ,n\ 

/3a k 

= 0 

Vz e [b2k-i,a 2 (n-k)+i] 

k e [l,n] 

Oik 

_ k 

A(l+2n) 

y Z € [a 2 (n-jfc)+l, a 2fc]) 

k e [1 ,n] 


(108) 


which simplifies to the solution found previously in [l] for the case of a vanishing FI- 
parameter. 


A.2 Case II: n - \ 4^1, - \ < n < - b. 

L 2m J ’ m 2 m 2 

In this case, the point b 2n +i lies inside the interval \-A,B], but the point a 2 n+i does not. 
The ordering of the interior resonance points in this case will be as in the previous, but 
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with odd- and even 6-resonances interchanged, that is: 

-A < b\ < b 2n < «2n-l < 

< a 2 < 63 < 6 2 ( n _ 1) < a 2 ( n -i)-i < ... < 

< a 2n _ 2 < 6 2n _i < b 2 < a± < a 2n < b 2n+ 1 < B. 


( 109 ) 


Again, the density will be given by a piecewise constant one, which, using the notations 
a/Sfc,/3fc,/Sccfc, are governed by the following equations: 

1 


+ (3cx n —] t + /3o: n+ i_fc 

fin—l—k fin-k ~ 


( y .( 3 n — & + ck/^ 


n+l-fc — ^ 

2 Oik + «n-fc + «n+1-fc = 
together with the boundary conditions 

«o = /3-i = /3a 0 = /3a n+ i = 0. 
This gives us the eigenvalue density as 


k e [0,n] 
k e [0,n] 
k e [1 ,n\ 
k e [l,n], 


p(z) = 


a Pk 2 \ 

P k = A(3+2n) 

/3a k = 0 

k 

ak = A(l+2 n) 


Vz € [a 2 fc,6 2 fc + i ] k e [0,n] 

Vz e \b' 2 k+iyb 2 ( n —k )] k e [0,n] 

VT e [6 2(n -fc+l)t a 2(n-fc)+l] ^ e [1,™] 

e [a 2 (n-fc) + i,a 2 fc], fce[l, n] 


( 110 ) 


( 111 ) 


( 112 ) 


A.3 Case III: rc= [4^1, 

L 2 to J ’ to 2 to 2 

In this case, both a 2n +i and & 2n+1 will lie in the interior of the interval, and the ordering 
amongst the interior resonance points will be as in the second case, but with interchanged 
ordering amongst the a^’s. Furthermore, since the point a 2n+ 1 now enters the interior of 
the interval, this will be the leftmost resonance point, and not b\ as in the previous case. 
The ordering may be seen in figure [7J 


Once more, we write down a piecewise constant Ansatz for the density, the only 
change from the previous case is the numbering and precise appearance on the boundary 
conditions for cc*,, (3k,oi(3k, P&k. Let our Ansatz be: 


Otk 

Vz e [a 2 fc,a 2 ( n _fc) + i], 

k e [0,n] 

a/3 k 

yz e [a 2 ( n _A,) + i,6 2 fc + i] 

k e [0,n] 

Pk 

yz e [^2A;+1^2(n-fc)] 

k e [0,77,] 

fiotk 

yz e [6 2 ( n+ i_/ c ),a 2 / c ] 

k e [ 1 , 77 ] 
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( 113 ) 











Figure 7: Resonance points originating from the interval endpoint —A and B (which are shown 
in blue). Resonances originating from —A are shown in magenta (even, a 2 fc), and purple (odd 
resonances, a 2 fc_i). Similarly, resonances originating from the rightmost interval endpoint, B, 
are shown in yellow (even, 6 2 fc) and red (odd, & 2 fc+i)- These resonances exist for all k < 4^21. 
Furthermore, the resonances under the maps 2 ->■ —z — m\ and z -*■ m 2 — z are shown underneath. 


Then, the saddle-point equation (18) takes the form: 


2 Oik "T OLn-k ■*" ®re+l —k 

2tr/3fc fiotn-k "t floin+l-k — 

2 /3k + fin—l—k Pn-k — 


2 f3a k + aj3 n _ k + a/3, 


n+l-k 


A 


k € [0,77.] 
k <E [0,77] 
k e [0,?7] 
k e [ 1 , 77 ], 


with 


p(z) = 


(114) 


a 2n+ i = P-i = Pa 0 = /3a n+1 = 0. (115) 

Hence, in general, the eigenvalue density in this case is given by: 

= x(3+2n) 6 [a2fc,a 2 (n-fc)+l], k 6 [0,77] 

apk = Jx e [ a 2 (n-k)+i,hk+i] k e [0, n\ 

Pk = A(3+2n) 6 [hk+i,h(n-k)] A; e [0,77] 

floik = 0 \/z € [b 2 (n+i-k),a 2 k\ k e [l,n] 

which indeed reduces to the solution previously found for vanishing Fl-parameter as B 
A and ( -*■ 0. 


(116) 
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A.4 Determining the interval endpoints 


Once more, the rightmost interval endpoint may be determined from the integral equation 
(16), where the two first integrals are simply determined from the normalisation condition. 
However, the third one is not. The sign-function in that integral will take the value one 
for y > m 2 - B = b\, whereas it will take the value -1 for y < m 2 - B = 61 , and so we find 
B as: 


B 


(3 - J^ Ayp(y) * J b 1 iyp(y))- 


(117) 


These integrals will depend on the ordering on the resonance points and thus will take 
different values in the different cases. Define the distances between one a and one b- 
resonance point as ab, the distance between two a-resonances as a, and between two 
6 -resonances as b. The integral may then be divided as: 


B =- ( 3 - ab Y oifik - a X! - b X! fa 

v <61 <61 <61 

+ ab Y, a fa + a Y a k + *>£&)■ 


(118) 


>6l >6l >6l 

Inserting the expressions for the eigenvalue densities, one then finds: 

£? =! ( 3 + ab 11 

2 V 


2A 


(119) 


-aj>*-b£/3 fc + a£ OLk + 

<6l <6l >6l >6l 


Furthermore, the normalisation condition on the eigenvalue density will give us another 
condition on the interval endpoints. In the context of ab,a and b, this takes the form: 


ab Y, a Pk + a Y a k + b X! Pk =1 


<=> 

ab ^ 7 - 7 - + a ^ ay. + b ^ = 1 . 

k k 


( 120 ) 


2A 


Together, the equations (119) and (120) gives us two relations between the interval end¬ 
points and the parameters m, ( and A, and these may be used to determine the interval 
endpoints in terms of these quantities. 

These expressions depend on the precise appearance of the eigenvalue densities, as 
well as on the ordering of the resonance points around b\. Therefore they will give rise 
to different results for the interval endpoints in the three cases. Inserting the appropriate 
limits in the sum, and the corresponding number of points to the left/right of b\, one finds 
the interval endpoints described in (47) - (49). 
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